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1.  Introduction 


Croconic  acid  (CA,  C2H2O5,  4,5-dihydroxycyclopentenetrione)  is  a  member  of  the 
family  of  oxycarbon  acids,  which  includes  deltic  acid  (C3H2O2),  squaric  acid 
(C4H2O4),  and  rhodizonic  acid  (C6H2O6).  The  molecular  structure  of  CA  is  shown 
in  Fig.  1.  Recent  research  conducted  at  the  US  Army  Research  Laboratory 
identified  a  high-pressure  polymorph  of  CA  (P-CA),  which  has  been  suggested  to 
have  high-energy  density  and  may  have  applications  as  a  propellant.1 


HO  OH 


O 


Fig.  1  Molecular  structure  of  CA  (This  image  is  copied  from  Wikipedia  in  the  public 
domain  and  is  free  of  copyright  and  can  be  used  without  any  license  requirements. 
https://en.wikipedia.org/wiki/Croconic_acid.) 

Large-scale  production  of  P-CA  is  hindered  by  the  high  pressures  necessary  to 
induce  the  a-  to  P-phase  change.  In  a  Paris  Edinburgh  large-volume  press,  the 
p-CA  polymorph  was  synthesized  through  the  compression  of  as-received,  a-CA 
to  9  GPa  for  a  minimum  of  24  h;  slightly  lower  transition  pressures  have  been  noted 
in  the  diamond  anvil  cell.2  The  necessity  of  the  high  pressure  in  its  synthesis  limits 
the  annual  production  of  p-CA,  using  one  Paris  Edinburgh  cell,  to  approximately 
13  g/yr.1  Previous  research  investigated  a  potential  means  of  increasing  production 
via  water-recrystallization.1  Such  methods  were  found  to  be  unsuccessful,  with  the 
initial  color  shift  of  CA,  which  was  previously  determined  to  be  indicative  of  the 
a-  to  P-phase  change,  found  to  be  a  particle-size  effect  and  not  a  result  of  the  desired 
phase  change. 

This  report  investigates  efforts  to  synthesize  P-CA  through  ball  milling.  In  ball 
milling,  high  pressures  are  generated  by  the  impact  of  a  ball  with  a  vessel  wall  while 
the  vessel  is  agitated.  This  method  is  well  established  in  the  literature  as  a  technique 
to  induce  mechanical  alloying  in  the  production  of  nano-crystalline  materials3  and 
reactive  materials.4  It  was  hoped  that  the  dynamic  pressures  generated  during 
milling  would  allow  for  rapid  production  of  P-CA,  but  this  was  found  to  not  be  the 
case. 
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2.  Experimental 


2.1  Ball  Milling 

Several  kinds  of  ball  mills  are  commonly  used  in  material  processing  including 
attritor  mills,  rolling  ball  mills,  and  planetary  ball  mills.  This  study  used  a  vibratory 
ball  mill;  such  mills  consist  of  a  vessel  containing  a  ball  that  is  moved  rapidly  back 
and  forth,  resulting  in  the  ball  colliding  with  the  ends  of  the  vessel  on  each  swing. 

A  Wig-L-Bug  (Dentsply  Rinn)  ball  mill  was  used  to  induce  high  dynamic  pressures 
in  as-received  CA  powder.  The  Wig-L-Bug  is  similar  in  many  respects  to  the 
commonly  used  “SPEX-type”  high-energy  ball  mill  but  on  a  smaller  scale.  The 
vessel  and  ball  in  the  Wig-L-bug  are  both  stainless  steel.  The  ball  has  a  diameter  of 
6.35  mm  and  a  mass  of  1  g;  the  vessel  outer  diameter  is  12.7  mm,  and  the  vessel 
length  is  27.7  mm.  A  series  of  milling  experiments  were  conducted,  varying  the 
ball  to  sample  mass  ratio  and  the  amount  of  milling  time.  Experiments  conducted 
are  summarized  in  Table  1. 


Table  1  Parameters  for  milling  experiments 


Sample  no. 

Ball-to-sample  mass 
ratio 

Mill  time 
(min) 

0  (control) 

NA 

0 

1 

20:1 

1 

2 

10:1 

1 

3 

5:1 

1 

4 

10:1 

5 

5 

10:1 

10 

NA  =  not  applicable 


2.2  High-Speed  Video 

A  Fastcam  SA-Z  (Photron)  high-speed  digital  video  camera  was  used  to  record  the 
motion  of  the  Wig-L-Bug  ball  mill;  from  this  motion  it  was  possible  to  determine 
the  velocity  parameters  needed  for  Hertzian  stress  calculations.  A  small  piece  of 
tape  was  placed  on  the  mill  to  aid  in  motion  tracking.  The  mill  was  activated  and 
then  the  camera  was  triggered  to  record  its  motion.  The  frame  rate  of  the  camera 
was  5000  fps  and  the  exposure  time  was  199  ps. 
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2.3  Differential  Scanning  Calorimetry 


Differential  scanning  calorimetry  (DSC)  experiments  were  conducted  to  determine 
how  milling  effected  the  CA  powder.  A  Q2000  DSC  (TA  Instruments)  was  used  to 
heat  a  small  sample  of  material  at  10  °C  /min  from  ambient  temperature  to  475  °C. 
During  the  heating  process,  the  furnace  was  purged  with  nitrogen  flowing  at 
100  mL/min.  Three  DSC  traces  were  collected  for  each  sample. 

3.  Results  and  Discussion 


3.1  Milling  Pressures 

The  high  pressures  achievable  in  this  type  of  mill  result  from  stresses  that  develop 
in  the  milled  material  as  it  is  trapped  between  the  ball  and  end  wall.  These  contact 
stresses  were  first  studied  by  Hertz5  and  referred  to  as  Hertzian  stresses.6  Maurice 
and  Courtney7  applied  the  concept  of  Hertzian  stresses  to  high-energy  ball  milling 
to  determine  the  parameters  important  to  mechanical  alloying.  They  determined  an 
expression  for  the  maximum  pressure  achieved  in  a  mill,  given  in  Eq.  1  as 


(1) 


and  an  expression  for  the  total  duration  of  the  ball-wall  impact  event  (2 r),  given  in 
Eq.  2  as 


2t  —  gTv  °-2{p/Eeff)°AR  , 


(2) 


where  pmax  is  the  maximum  pressure,  gP  and  g T  are  geometric  constants,  v  is  the 
relative  ball-wall  velocity,  p  is  the  ball  density,  Eeff  is  the  effective  tensile  modulus 
of  the  ball  and  container,  and  R  is  the  ball  radius.  In  this  study,  the  ball  and  container 
are  both  made  of  the  same  material,  so  the  effective  modulus  is  the  tensile  modulus 
of  stainless  steel.  The  values  of  these  parameters  are  listed  in  Table  2.  The  values 
of  gP  and  g T  are  for  a  ball  impacting  a  flat  surface,  which  was  deemed  to  be  an 
acceptable  approximation  given  the  large  difference  in  radius  of  curvature  between 
the  ball  and  vessel  ends. 
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Table  2  Values  for  Eqs.  1  and  27-x 


Parameter 

Value 

Unit 

9p 

0.3521 

9 T 

6.4034 

P 

7450 

kg/m3 

Eeff 

210  x  109 

Pa 

R 

3.18  x  103 

m 

To  determine  the  impact  velocity  of  the  ball  on  the  end  wall  of  the 
Wig-L-Bug,  a  high-speed  video  of  the  mill  in  motion  was  imported  into  ImageJ 
(National  Institutes  of  Health)  for  image  analysis.  The  video  was  advanced  frame- 
by-frame  and  the  leading  edge  of  the  small  piece  of  tape  placed  on  the  mill  clamp 
was  determined  in  each  frame.  A  screen  capture  from  ImageJ  after  motion  analysis 
is  shown  in  Fig.  2. 


Fig.  2  ImageJ  motion  analysis  of  Wig-L-Bug  during  operation.  White  crosses  represent  the 
position  of  the  center  of  the  Wig-L-Bug  clamp,  and  yellow  markings  denote  the  video  frame 
at  which  the  positon  was  determined. 
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The  coordinates  determined  by  motion  analysis  are  plotted  in  Fig.  3.  The  motion  of 
the  mill  is  3-dimensional  (with  the  z-component  not  recorded),  but  the  x-  and 
z-components  of  the  motion  are  significantly  smaller  than  the  y-component.  To 
simplify  analysis,  only  the  y-component  of  the  motion  was  used  to  determine  the 
relative  impact  velocity. 


Fig.  3  Plotted  coordinates  of  Wig-L-Bug  motion  from  Fig.  2.  The  y-movement  is 
significantly  greater  in  amplitude  than  the  x-movement.  The  z-movement  is  not  recorded  in 
this  video,  but  it  is  on  the  order  of  the  x-movement. 

The  relative  impact  velocity  was  calculated  using  the  relationship  developed  by 
Maurice  and  Courtney  for  a  SPEX-type  mill,  given  in  Eq.  3  as7 

v  —  2  x/t  (3) 

where  v  is  the  relative  impact  velocity,  x  is  the  distance  traveled  by  the  mill  in  one 
full  vibration  cycle  (i.e.,  twice  the  amplitude),  and  t  is  the  time  of  one  milling  cycle. 

Plotting  the  y-component  of  the  motion  of  the  mill  (Fig.  4),  the  frequency  of  the 
mill  is  approximately  50  Hz,  and  the  amplitude  of  the  motion  is  approximately 
2  cm.  Thus  x  =  0.04  m  and  t  =  0.02  s,  roughly.  When  Eq.  3  is  used,  this  results  in  a 
relative  impact  velocity  of  v  ~  4  m/s.  At  this  impact  velocity  the  peak  impact 
pressure  is  4.2  GPa,  and  the  impact  duration  is  1.6  qs. 
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Fig.  4  The  y-component  of  Wig-L-Bug  motion  in  time.  The  frequency  is  approximately 
50  Hz,  and  the  amplitude  is  approximately  2  cm. 

The  maximum  pressure  and  duration  of  impact  over  a  range  of  impact  velocities, 
calculated  using  Eqs.  1  and  2  with  the  parameters  listed  in  Table  2,  are  plotted  in 
Fig.  5.  The  velocity  needed  to  reach  9-GPa  peak  pressure  is  approximately 
28  m/s,  which  has  an  impact  duration  of  1.1  ps. 


Fig.  5  Peak  pressure  ( pmax )  and  duration  of  impact  (2t)  as  a  function  of  relative  impact 
velocity  for  a  stainless  steel  ball  and  vessel 
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3.2  Differential  Scanning  Calorimetry 


Representative  heat  flow  traces  showing  the  exothermic  decomposition  of  CA 
samples  are  shown  in  Fig.  6.  The  integrated  exotherm  areas  for  all  samples  are 
shown  in  Fig.  7.  From  these  figures,  it  is  apparent  that  varying  the  ball  to  sample 
mass  ratio  (Samples  1-3)  did  not  have  an  effect  on  the  resulting  exothermicity  of 
the  milled  material  compared  to  the  control  sample  (Sample  0),  and  extending  the 
milling  time  had  only  a  moderate  effect.  Samples  milled  for  5  min  (Sample  4) 
showed  somewhat  higher  exothermicity  than  the  control  sample  but  with  significant 
variation  in  the  results.  Samples  milled  for  10  min  (Sample  5)  showed  a  significant 
reduction  in  exothermicity,  again  with  large  variation.  This  may  be  due  to  sample 
degradation  caused  by  elevated  temperatures  induced  by  frictional  heating  in  the 
mill. 


Fig.  6  Representative  heat  flow  traces  for  milled  CA  samples.  Curves  have  been  offset  to 
enhance  clarity.  Discontinuities  in  the  traces  are  attributed  to  slight  movement  of  the  sample 
pans  due  to  gas  production  by  the  sample  and  do  not  significantly  affect  energy  release 
calculations. 
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Fig.  7  Energy  released  by  the  decomposition  of  CA,  determined  by  integrating  the  DSC 
trace.  Error  bars  represent  95%  confidence  intervals. 

4.  Conclusions 


From  the  DSC  measurements,  it  is  apparent  that  the  currently  available  milling 
equipment  is  likely  not  an  effective  method  to  induce  the  desired  phase  change  ot¬ 
to  fl-CA.  This  agrees  with  the  pressures  and  impact  durations  calculated  based  on 
high-speed  video  analysis  of  the  Wig-L-Bug.  The  pressures  achievable  with  this 
small  mill  are  less  than  half  that  required  to  induce  the  phase  change.  These 
pressures  are  also  only  sustained  for  times  on  the  order  of  microseconds,  rather  than 
the  several  hours  long  duration  currently  employed. 

Higher  pressures  could  be  achieved  by  increasing  the  relative  impact  velocity  of 
the  mill  or  by  modifying  the  mill/ball  materials.  However,  due  to  the  very  short 
impact  durations  of  ball  milling  it  is  likely  that  pressures  significantly  greater  than 
9  GPa  would  be  required.  It  will  also  be  important  to  minimize  heat  generation 
caused  by  milling  to  prevent  material  decomposition.  This  may  be  possible  through 
the  use  of  cryogenic  milling. 


Sample  0  Sample  1  Sample  2  Sample  3  Sample  4  Sample  5 
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